Genomic and biological characterization of bacteriophages against Enterobacter cloacae, a high-priority pathogen.

Enterobacter cloacae is a clinically significant pathogen due to its multi-resistance to antibiotics, presenting a challenge in the treatment of infections. As concerns over antibiotic resistance escalate, novel therapeutic approaches have been explored. Bacteriophages, characterized by their remarkable specificity and ability to self-replicate within target bacteria, are emerging as a promising alternative therapy. In this study, we isolated and partially characterized nine lytic bacteriophages targeting E. cloacae, with two selected for comprehensive genomic analysis based on their host range and bacteriolytic activity. All identified phages exhibited a narrow host range, demonstrated stability within a temperature range of 30-60 °C, displayed pH tolerance from 3 to 10, and showed an excellent bacteriolytic capacity for up to 18 h. Notably, the fully characterized phage genomes revealed an absence of lysogenic, virulence, or antibiotic-resistance genes, positioning them as promising candidates for therapeutic intervention against E. cloacae-related diseases. Nonetheless, translating this knowledge into practical therapeutic applications mandates a deeper understanding of bacteriophage interactions within complex biological environments.

Phylogenomic analysis of the Ponticus clade: strains isolated from the spotted rose snapper (Lutjanus guttatus).

The family Vibrionaceae is classified into many clades based on their phylogenetic relationships. The Ponticus clade is one of its clades and consists of four species, Vibrio panuliri, V. ponticus, V. rhodolitus, and V. taketomensis. Two strains, CAIM 703 and CAIM 1902, were isolated from the diseased spotted rose snapper external lesion (Lutjanus guttatus), they were analyzed to determine their taxonomic position, a phylogenetic analysis was performed based on the 16S rRNA sequences proved that the two strains are members of the genus Vibrio and they belong to the Ponticus clade. Then, a phylogenomic analysis was performed with four type strains and four reference strains isolated from marine organisms and aquatic environments. Multilocus Sequence Analysis (MLSA) of 139 single-copy genes showed that CAIM 703 and CAIM 1902 belong to V. panuliri. The 16S rRNA sequence similarity value between CAIM 703 and CAIM 1902 was 99.61%. The Ponticus clade species showed Average Nucleotide Identity (ANI) values between 78 to 80% against the two strains for ANIb, except V. panuliri LBS2T (99% and 100% similarity). Finally, this analysis represents the first phylogenomic analysis of the Ponticus clade where V. panuliri strains are reported from Mexico.

Exploration of low-frequency allelic variants of SARS-CoV-2 genomes reveals coinfections in Mexico occurred during periods of VOCs turnover.

A total of 14 973 alleles in 29 661 sequenced samples collected between March 2021 and January 2023 by the Mexican Consortium for Genomic Surveillance (CoViGen-Mex) and collaborators were used to construct a thorough map of mutations of the Mexican SARS-CoV-2 genomic landscape containing Intra-Patient Minor Allelic Variants (IPMAVs), which are low-frequency alleles not ordinarily present in a genomic consensus sequence. This additional information proved critical in identifying putative coinfecting variants included alongside the most common variants, B.1.1.222, B.1.1.519, and variants of concern (VOCs) Alpha, Gamma, Delta, and Omicron. A total of 379 coinfection events were recorded in the dataset (a rate of 1.28 %), resulting in the first such catalogue in Mexico. The most common putative coinfections occurred during the spread of Delta or after the introduction of Omicron BA.2 and its descendants. Coinfections occurred constantly during periods of variant turnover when more than one variant shared the same niche and high infection rate was observed, which was dependent on the local variants and time. Coinfections might occur at a higher frequency than customarily reported, but they are often ignored as only the consensus sequence is reported for lineage identification.

Gestational diabetes mellitus: genetic factors, epigenetic alterations, and microbial composition.

Gestational diabetes mellitus (GDM) is a common metabolic disorder, usually diagnosed during the third trimester of pregnancy that usually disappears after delivery. In GDM, the excess of glucose, fatty acids, and amino acids results in foetuses large for gestational age. Hyperglycaemia and insulin resistance accelerate the metabolism, raising the oxygen demand, and creating chronic hypoxia and inflammation. Women who experienced GDM and their offspring are at risk of developing type-2 diabetes, obesity, and other metabolic or cardiovascular conditions later in life. Genetic factors may predispose the development of GDM; however, they do not account for all GDM cases; lifestyle and diet also play important roles in GDM development by modulating epigenetic signatures and the body's microbial composition; therefore, this is a condition with a complex, multifactorial aetiology. In this context, we revised published reports describing GDM-associated single-nucleotide polymorphisms (SNPs), DNA methylation and microRNA expression in different tissues (such as placenta, umbilical cord, adipose tissue, and peripheral blood), and microbial composition in the gut, oral cavity, and vagina from pregnant women with GDM, as well as the bacterial composition of the offspring. Altogether, these reports indicate that a number of SNPs are associated to GDM phenotypes and may predispose the development of the disease. However, extrinsic factors (lifestyle, nutrition) modulate, through epigenetic mechanisms, the risk of developing the disease, and some association exists between the microbial composition with GDM in an organ-specific manner. Genes, epigenetic signatures, and microbiota could be transferred to the offspring, increasing the possibility of developing chronic degenerative conditions through postnatal life.

SARS-CoV-2 Omicron variants BA.4 and BA.5 dominated the fifth COVID-19 epidemiological wave in Mexico.

In Mexico, the BA.4 and BA.5 Omicron variants dominated the fifth epidemic wave (summer 2022), superseding BA.2, which had circulated during the inter-wave period. The present study uses genome sequencing and statistical and phylogenetic analyses to examine these variants' abundance, distribution, and genetic diversity in Mexico from April to August 2022. Over 35 % of the sequenced genomes in this period corresponded to the BA.2 variant, 8 % to the BA.4 and 56 % to the BA.5 variant. Multiple subvariants were identified, but the most abundant, BA.2.9, BA.2.12.1, BA.5.1, BA.5.2, BA.5.2.1 and BA.4.1, circulated across the entire country, not forming geographical clusters. Contrastingly, other subvariants exhibited a geographically restricted distribution, most notably in the Southeast region, which showed a distinct subvariant dynamic. This study supports previous results showing that this region may be a significant entry point and contributed to introducing and evolving novel variants in Mexico. Furthermore, a differential distribution was observed for certain subvariants among specific States through time, which may have contributed to the overall increased diversity observed during this wave compared to the previous ones. This study highlights the importance of sustaining genomic surveillance to identify novel variants that may impact public health.

Genomic Surveillance of SARS-CoV-2 in México: Three Years since Wuhan, China's First Reported Case.

OBJECTIVE: The aim of this work was to analyze the metadata of the SARS-CoV-2 sequences obtained from samples collected in Mexico from 2020 to 2022. MATERIALS AND METHODS: Metadata of SARS-CoV-2 sequences from samples collected in Mexico up to 31 December 2022 was retrieved from GISAID and manually cured for interpretation. RESULTS: As of December 2022, Mexican health authorities and the scientific community have sequenced up to 81,983 SARS-CoV-2 viral genomes deposited in GISAID, representing 1.1% of confirmed cases. The number of sequences obtained per state corresponded to the gross domestic product (GDP) of each state for the first (Mexico City) and the last (Tlaxcala). Approximately 25% of the sequences were obtained from CoViGen-Mex, an interdisciplinary initiative of health and scientific institutions to collect and sequence samples nationwide. The metadata showed a clear dominance of sequences retrieved by women. A similar variant distribution over time was found in Mexico and overseas, with the Omicron variant predominating. Finally, the age group with the highest representation in the sequences was adults aged 21 to 50 years, accounting for more than 50% of the total. CONCLUSIONS: Mexico presents diverse sociodemographic and economic characteristics. The COVID-19 pandemic has been and continues to be a challenge for collaboration across the country and around the world.

First report of Pectobacterium polaris causing aerial stem rot of potato in Mexico.

In Mexico, potato (Solanum tuberosum L.) is one of the most important vegetable crops for local consumption and industry. More than 1.8 million tons of potatoes are produced annually, of which the state of Sinaloa contributes with 21.5% (SIAP. 2022). In January 2020, potato plants (cv. FL1867) showing aerial stem rot symptoms were observed in a commercial field from the Santa Rosa Valley, in Northern Sinaloa with an incidence of 36%. Putative pectolytic bacteria showing pitting on crystal violet pectate (CVP) plates were restreaked and purified onto Nutritive Agar (NA) medium at 28°C. Four independent isolates were obtained (L25F-83, L25F-105, L25F-115, and L25F-125) from four symptomatic stems with biochemical and morphological characteristics related to Pectobacterium, such as catalase positive, oxidase negative, pectinolytic activity, Gram-negative and non-fluorescent in B-King medium. Bacterial gDNA was used for amplification and sequencing of two housekeeping genes (dnaX and leuS) (Portier et al. 2019). The nucleotide sequence identity between our isolates was 100% with both housekeeping genes (dnaX, OP376536-OP376539 and leuS, OP376540-OP376543). The BLASTn analysis of dnaX gene shared 98.98% and 99.19% identity with two soft-rot-causing bacterial strains NIBIO1006T (CP017481) and NIBIO1392 (CP017482) of Pectobacterium polaris, respectively; and with leuS gene shared 99.56% identity with P. polaris type strain NIBIO1006T. To further validate the identification, two strains, S5 (isolate L25F-105) and S6 (L25F-125), were selected for whole genome sequencing (WGS). The ANI values for closely related species were calculated using the Orthologous Average Nucleotide Identity (Ortho-ANI) Software Tool (OAT) (Lee et al. 2016). The Type (Strain) Genome Server (TYGS) was used for accurate genome-based taxonomy (https://tygs.dsmz.de) (Meier-Kolthoff et al. 2019). The genomes of P. polaris strains S5 (4811345 pb, GC=52%, AULSZ000000000) and S6 (4809754 pb, GC=52%, JAULTA000000000) revealed 96.86% and 96.07% Ortho-ANI and 73.6% and 66.5% dDDH with P. polaris type strain NIBIO1006T and P. parvum strain CFBP8630, respectively. The MLSA was performed on concatenated complete sequences of dnaX (OR470476, OR470477), leuS (OR470484, OR470485), recA (OR470488, OR470488), acnA (OR470474, OR470475), gapA (OR470478, OR470479), gyrA (OR470480, OR470481), icdA (OR470482, OR470483), proA (OR470486, OR470487), and rpoA genes (OR470490, OR470491). The consensus tree, constructed using the maximum likelihood method (MEGA 7.0), clustered strains S5 and S6 with P. polaris strains NIBIO1006T and NIBIO1392. The four isolates resulted pathogenic in tuber slices and potato seedlings (cv. Fianna) 24 and 72 h, respectively, after being inoculated with 30 µL bacterial suspension (1 X 108 CFU/ml) and incubated at 28 °C and 85% relative humidity. Bacterial colonies were reisolated from the affected tissue and identified with the same PCR primers as described above. Accordingly, P. polaris isolates S5 and S6, fulfill Koch's postulates for aerial stem rot of potato. To our knowledge, this is the first report of P. polaris causing aerial stem rot of potato in Mexico. This bacterium could be a significant threat to the local potato producers; therefore, an accurate and sensitive method of detection and epidemiological studies are needed to support an effective disease diagnosis and management program.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico.

Over 200 different SARS-CoV-2 lineages have been observed in Mexico by November 2021. To investigate lineage replacement dynamics, we applied a phylodynamic approach and explored the evolutionary trajectories of five dominant lineages that circulated during the first year of local transmission. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in Mexico. Lineages B.1.1.222 and B.1.1.519 exhibited similar dynamics, constituting clades that likely originated in Mexico and persisted for >12 months. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. For the largest B.1.617.2 clades, we further explored viral lineage movements across Mexico. Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

Vibrio eleionomae sp. nov., isolated from shrimp (Penaeus vannamei) pond water.

A novel Vibrio strain (CAIM 722T=SW9T=DSM 24596T) was isolated in 2003 from water of a shrimp (Penaeus vannamei) culture pond located in Los Mochis, Sinaloa, Mexico, and taxonomically characterized using a polyphasic approach. The 16S rRNA gene sequence clustered within those of the genus Vibrio, showing high similarity to the type strains of the Porteresiae clade. Multilocus sequence analysis using eight housekeeping genes (ftsZ, gapA, gyrB, mreB, pyrH, recA, rpoA, topA and 16S rRNA) and phylogenetic analysis with 139 single-copy genes showed that the strain forms an independent branch. Whole genome sequencing and genomic analyses (average nucleotide identity, OrthoANI, average amino acid identity and in silico DNA-DNA hybridization) produced values well below the thresholds for species delineation with all methods tested. In addition, a phenotypic characterization was performed to support the description and differentiation of the novel strain from related taxa. The results obtained demonstrate that the strain represent a novel species, for which the name Vibrio eleionomae sp. nov. is proposed.

SARS-CoV-2 BW lineage, a fast-growing Omicron variant from southeast Mexico bearing relevant escape mutations.

PURPOSE: The swift expansion of the BW.1 SARS-CoV-2 variant coincided with a rapid increase of COVID-19 cases occurring in Southeast Mexico in October, 2022, which marked the start of Mexico's sixth epidemiological wave. In Yucatan, up to 92% (58 of 73) of weekly sequenced genomes between epidemiological week 42 and 47 were identified as either BW.1 or its descendant, BW.1.1 in the region, during the last trimester of 2022. In the current study, a comprehensive genomic comparison was carried out to characterize the evolutionary history of the BW lineage, identifying its origins and its most important mutations. METHODS: An alignment of all the genomes of the BW lineage and its parental BA.5.6.2 variant was carried out to identify their mutations. A phylogenetic and ancestral sequence reconstruction analysis with geographical inference, as well as a longitudinal analysis of point mutations, were performed to trace back their origin and contrast them with key RBD mutations in variant BQ.1, one of the fastest-growing lineages to date. RESULTS: Our ancestral reconstruction analysis portrayed Mexico as the most probable origin of the BW.1 and BW.1.1 variants. Two synonymous substitutions, T7666C and C14599T, support their Mexican origin, whereas other two mutations are specific to BW.1: S:N460K and ORF1a:V627I. Two additional substitutions and a deletion are found in its descending subvariant, BW.1.1. Mutations found in the receptor binding domain, S:K444T, S:L452R, S:N460K, and S:F486V in BW.1 have been reported to be relevant for immune escape and are also key mutations in the BQ.1 lineage. CONCLUSIONS: BW.1 appears to have arisen in the Yucatan Peninsula in Southeast Mexico sometime around July 2022 during the fifth COVID-19 wave. Its rapid growth may be in part explained by the relevant escape mutations also found in BQ.1.

Omicron-BA.1 Dispersion Rates in Mexico Varied According to the Regional Epidemic Patterns and the Diversity of Local Delta Subvariants.

PURPOSE: The Omicron subvariant BA.1 of SARS-CoV-2 was first detected in November 2021 and quickly spread worldwide, displacing the Delta variant. In this work, a characterization of the spread of this variant in Mexico is presented. METHODS: The time to fixation of BA.1, the diversity of Delta sublineages, the population density, and the level of virus circulation during the inter-wave interval were determined to analyze differences in BA.1 spread. RESULTS: BA.1 began spreading during the first week of December 2021 and became dominant in the next three weeks, causing the fourth COVID-19 epidemiological surge in Mexico. Unlike previous variants, BA.1 did not exhibit a geographically distinct circulation pattern. However, a regional difference in the speed of the replacement of the Delta variant was observed. CONCLUSIONS: Viral diversity and the relative abundance of the virus in a particular area around the time of the introduction of a new lineage seem to have influenced the spread dynamics, in addition to population density. Nonetheless, if there is a significant difference in the fitness of the variants, or if the time allowed for the competition is sufficiently long, it seems the fitter virus will eventually become dominant, as observed in the eventual dominance of the BA.1.x variant in Mexico.

Characterization and genome analysis of six novel Vibrio parahaemolyticus phages associated with acute hepatopancreatic necrosis disease (AHPND).

Vibrio parahaemolyticus causes acute hepatopancreatic necrosis disease (AHPND) in farmed shrimp. Due to its damage potential, which could be as high as a 100% mortality rate, bacteriophages have emerged as a promising natural control intervention other than antibiotics, yet multiple roadblocks need to be overcome. In this study, six bacteriophages isolated from seafood samples, seawater, and estuary water in Sinaloa, Mexico, demonstrated a narrow host range among Mexican AHPND-causing V. parahaemolyticus. All bacteriophages are composed of a double-stranded DNA genome with lengths ranging between 43,268 and 57,805 bp. All six phages exhibited latency periods of 10-30 min and burst sizes of 34-168 viral particles per infected cell. The optimal MOI for bacteriophage propagation was 0.01-1. No transfer RNA (tRNA), virulence, or resistance genes were found in either genome, and the life cycle of these phages was classified as virulent by the PhageAI platform. Phylogenetic and comparative genomics analyzes assigned phages M3, C2, M9, and M83 as new species not yet reported within the genus Maculvirus, Autographiviridae family. ALK and CHI phages were assigned as new members of a new genus not yet classified within the subfamily Queuovirinae. The findings highlight the potential of CHI, ALK, M3, C2, M9, and M83 as promising alternatives against AHPND-causing V. parahaemolyticus from Mexico.

A Broad-Host-Range Phage Cocktail Selectively and Effectively Eliminates Vibrio Species from Shrimp Aquaculture Environment.

The protective effects of a phage cocktail composed of vB_Vc_SrVc2 and vB_Vc_SrVc9 were tested in Pacific white shrimp (Litopenaeus vannamei) postlarvae, which were originally isolated from diseased shrimps and selected due to their broad-host-range properties against several pathogenic Vibsrio species. We used culture-dependent and culture-independent approaches to explore its effect on bacterial communities associated with shrimp postlarvae. Both methods revealed that the levels of Vibrio species were significantly reduced after phage cocktail administration. Phage-treated shrimp also exhibisuppted lesser damage and higher lipid accumulation in B cells of the hepatopancreas, as revealed by histopathological examination. Taken together, this study provides clear evidence that phage therapy can selectively and effectively reduce Vibrio species, thereby providing an environmentally safe alternative to the prophylactic use of antibiotics in shrimp aquaculture.

Genomic Characterization of Twelve Lytic Bacteriophages Infecting Midgut Bacteria of Aedes aegypti.

Mosquito-borne diseases such as malaria and dengue are global severe public health threats. Due to the lack of efficient control methods, alternative approaches to decreasing arboviral transmitted diseases are prioritized to reduce morbidity and mortality in every endemic region. Mosquito midgut bacteria play an essential role in physiological development, fitness, and the arthropods´ vectorial capacity. Bacteriophages are viruses that infect bacteria and are considered a promising biocontrol method by eliminating midgut microbiota that plays an essential role in mosquitoes´ health. Here, we isolate and identify 22 bacteria from mosquito´s midgut belonging to the genera Mesobacillus, Enterobacter, Klebsiella, Microbacterium, Micrococcus, Pantoea, Serratia, and Staphylococcus, mainly. Twelve phages with lytic activity against Enterobacter, Klebsiella, and Pantoea were also isolated. All 12 phages showed a double-stranded DNA genome, ranging from 36,790 to 149,913 bp, and were taxonomically classified as members of the Drexlerviridae family, Molineuxvirinae, Studiervirinae, and Vequintavirinae subfamilies. Open reading frames associated with phage structure, packing, host lysis, DNA metabolism, and additional functions were predicted in all 12 phage genomes, while tRNAs were predicted in five phage genomes. In addition, the life cycle was predicted as virulent for the 12 phages, and no antibiotic resistance, virulence, allergenic, or lysogenic genes were found in either genome. These findings suggest that the 12 phages have biocontrol potentials; however, it is necessary to elucidate specific bacterial host's roles and then the phages' ability to serve as effective vector control.

Draft genome sequences of Bacillus pumilus 36R ATNSAL and B. safensis 13L LOBSAL, two potential candidate probiotic strains for shrimp aquaculture.

OBJECTIVES: This work aimed to isolate bacterial strains with antagonist activity against Vibrio parahaemolyticus, the causative agent of acute hepatopancreatic necrosis disease (VPAHPND) that was isolated from outbreaks in Mexico. Here, we report the draft genome sequences of two antagonistic strains, isolated from saline sediment in Sonora, Mexico. METHODS: Cross-streak and well diffusion tests were employed to find the bacterial strains with higher inhibitory activity against VPAHPND. The whole genomes of B. pumilus 36R ATNSAL and B. safensis 13L LOBSAL were sequenced using Ion TorrentTM (PGM) and Illumina MiseqTM platforms, respectively. Annotation was performed using the RAST server, and the genes involved in the biosynthesis of bacterial secondary metabolites were predicted using antiSMASH. RESULTS: Two bacterial isolates, B. safensis 13L LOBSAL and B. pumilus 36R ATNSAL, were chosen based on their strong antagonistic profiles. The genome of 36R ATNSAL was 3.94 Mbp in length and contained 3824 genes and a total of 4116 coding sequences (CDSs); the genome of 13L LOBSAL was 3.68 Mbp and contained 3619 genes and 3688 CDSs. Twenty-eight and 32 biosynthetic gene clusters responsible for putative antimicrobial metabolite production were identified in 36R ATNSAL and 13L LOBSAL, respectively. CONCLUSIONS: The two strains 13L LOBSAL and 36R ATNSAL showed excellent probiotic profiles in vitro. The genome sequences will help with the mining and reconstruction of metabolic pathways in Bacillus strains. Genome sequence-guided strain improvement could augment the probiotic potential of Bacillus strains for applications in shrimp aquaculture.

Global m6A RNA Methylation in SARS-CoV-2 Positive Nasopharyngeal Samples in a Mexican Population: A First Approximation Study.

The Severe Acute Respiratory Syndrome-Coronavirus-2 (SARS-CoV-2) is the causal agent of COVID-19 (Coronavirus Disease-19). Both mutation and/or recombination events in the SARS-CoV-2 genome have resulted in variants that differ in transmissibility and severity. Furthermore, RNA methylation of the N6 position of adenosine (m6A) is known to be altered in cells infected with SARS-CoV-2. However, it is not known whether this epitranscriptomic modification differs across individuals dependent on the presence of infection with distinct SARS-CoV-2 variants, the viral load, or the vaccination status. To address this issue, we selected RNAs (n = 60) from SARS-CoV-2 sequenced nasopharyngeal samples (n = 404) of 30- to 60-year-old outpatients or hospitalized individuals from the city of Mazatlán (Mexico) between February 2021 and March 2022. Control samples were non-infected individuals (n = 10). SARS-CoV-2 was determined with real-time PCR, viral variants were determined with sequencing, and global m6A levels were determined by using a competitive immunoassay method. We identified variants of concern (VOC; alpha, gamma, delta, omicron), the variant of interest (VOI; epsilon), and the lineage B.1.1.519. Global m6A methylation differed significantly across viral variants (p = 3.2 × 10-7). In particular, we found that m6A levels were significantly lower in the VOC delta- and omicron-positive individuals compared to non-infected individuals (p = 2.541236 × 10-2 and 1.134411 × 10-4, respectively). However, we uncovered no significant correlation between global m6A levels and viral nucleocapsid (N) gene expression or age. Furthermore, individuals with complete vaccination schemes showed significantly lower m6A levels than unvaccinated individuals (p = 2.6 × 10-4), and differences in methylation levels across variants in unvaccinated individuals were significant (p = 3.068 × 10-3). These preliminary results suggest that SARS-CoV-2 variants show differences in global m6A levels.

Prevalence and Genomic Diversity of Salmonella enterica Recovered from River Water in a Major Agricultural Region in Northwestern Mexico.

Salmonella enterica is a leading cause of human gastrointestinal disease worldwide. Given that Salmonella is persistent in aquatic environments, this study examined the prevalence, levels and genotypic diversity of Salmonella isolates recovered from major rivers in an important agricultural region in northwestern Mexico. During a 13-month period, a total of 143 river water samples were collected and subjected to size-exclusion ultrafiltration, followed by enrichment, and selective media for Salmonella isolation and quantitation. The recovered Salmonella isolates were examined by next-generation sequencing for genome characterization. Salmonella prevalence in river water was lower in the winter months (0.65 MPN/100 mL) and significantly higher in the summer months (13.98 MPN/100 mL), and a Poisson regression model indicated a negative effect of pH and salinity and a positive effect of river water temperature (p = 0.00) on Salmonella levels. Molecular subtyping revealed Oranienburg, Anatum and Saintpaul were the most predominant Salmonella serovars. Single nucleotide polymorphism (SNP)-based phylogeny revealed that the detected 27 distinct serovars from river water clustered in two major clades. Multiple nonsynonymous SNPs were detected in stiA, sivH, and ratA, genes required for Salmonella fitness and survival, and these findings identified relevant markers to potentially develop improved methods for characterizing this pathogen.

Dominance of Three Sublineages of the SARS-CoV-2 Delta Variant in Mexico.

In this study, we analyzed the sequences of SARS-CoV-2 isolates of the Delta variant in Mexico, which has completely replaced other previously circulating variants in the country due to its transmission advantage. Among all the Delta sublineages that were detected, 81.5 % were classified as AY.20, AY.26, and AY.100. According to publicly available data, these only reached a world prevalence of less than 1%, suggesting a possible Mexican origin. The signature mutations of these sublineages are described herein, and phylogenetic analyses and haplotype networks are used to track their spread across the country. Other frequently detected sublineages include AY.3, AY.62, AY.103, and AY.113. Over time, the main sublineages showed different geographical distributions, with AY.20 predominant in Central Mexico, AY.26 in the North, and AY.100 in the Northwest and South/Southeast. This work describes the circulation, from May to November 2021, of the primary sublineages of the Delta variant associated with the third wave of the COVID-19 pandemic in Mexico and highlights the importance of SARS-CoV-2 genomic surveillance for the timely identification of emerging variants that may impact public health.

The Alpha Variant (B.1.1.7) of SARS-CoV-2 Failed to Become Dominant in Mexico.

During the coronavirus disease 2019 (COVID-19) pandemic, the emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in the United Kingdom in September 2020, was well documented in different areas of the world and became a global public health concern because of its increased transmissibility. The B.1.1.7 lineage was first detected in Mexico during December 2020, showing a slow progressive increase in its circulation frequency, which reached its maximum in May 2021 but never became predominant. In this work, we analyzed the patterns of diversity and distribution of this lineage in Mexico using phylogenetic and haplotype network analyses. Despite the reported increase in transmissibility of the B.1.1.7 lineage, in most Mexican states, it did not displace cocirculating lineages, such as B.1.1.519, which dominated the country from February to May 2021. Our results show that the states with the highest prevalence of B.1.1.7 were those at the Mexico-U.S. border. An apparent pattern of dispersion of this lineage from the northern states of Mexico toward the center or the southeast was observed in the largest transmission chains, indicating possible independent introduction events from the United States. However, other entry points cannot be excluded, as shown by multiple introduction events. Local transmission led to a few successful haplotypes with a localized distribution and specific mutations indicating sustained community transmission. IMPORTANCE The emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) throughout the world were due to its increased transmissibility. However, it did not displace cocirculating lineages in most of Mexico, particularly B.1.1.519, which dominated the country from February to May 2021. In this work, we analyzed the distribution of B.1.1.7 in Mexico using phylogenetic and haplotype network analyses. Our results show that the states with the highest prevalence of B.1.1.7 (around 30%) were those at the Mexico-U.S. border, which also exhibited the highest lineage diversity, indicating possible introduction events from the United States. Also, several haplotypes were identified with a localized distribution and specific mutations, indicating that sustained community transmission occurred in the country.